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SUMMARY 

Bifurcation theory is used to analyse the nonlinear dynamic stability characteristics of single-degree- 
of-freedom motions of an aircraft or a flap about a trim position. The bifurcation theory analysis reveals 
that when the bifurcation parameter, e.g., the angle of attack, is increased beyond a orltical value at 
which the aerodynamic damping vanishes, a new solution representing finite-amplitude periodic motion 
bifurcates from the previously stable steady motion. The sign of a simple criterion, cast in terms of 
aerodynamic properties, determines whether the bifuroating solution is stable (supercritical) or unstable 
(subcritical) . For the pitching motion of a flat-plate airfoil flying at supersonio/hypersonic speed, and 
for oscillation of a flap at transonic speed, the bifurcation is subcritical, implying either that exchangee 
of stability between steady and periodic motion are accompanied by hysteresis phenomena, or that potentially 
large aperiodic departures from steady motion may develop. On the other hand, for the rolling oscillation 
of a slender delta wing in subsonic flight (wing rook), the bifurcation is found to be supercritical. This 
a«.d the predicted amplitude of the bifurcation periodic motion are in good agreement w. ¥ .t experiments. 


1 . INTRODUCTION 

Problems of aerodynamic stability of aircraft flying at small angles of attack have been studied exten- 
sively. With increasing angles of attack the problems become more complicated and typically involve non- 
linear phenomena such as coupling between modes, amplitude and frequency effects, and hysteresis. The need 
for investigating stability characteristics at high angles of attack was clearly demonstrated by Orlik- 
RUckemann (Ref. i) in his survey paper which deals largely with experiments. 

On the theoretical side, much of an extensive body of work is based on the linear theory, in which the 
unsteady flow is regarded as a small perturbation of some known steady flow (possibly nonlinear in, e.g., 
the angle of attack) that prevails under certain flight conditions. The question of the validity and limi- 
tations of such a linearized perturbation theory is of fundamental importance. Yet, it has rarely been 
investigated. One may argue that, in principle, it should be possible to advance to higher and higher 
angles of attack o by a series of linear perturbations since the solution at each step should Include a 
steady-state part which, when added to the previous steady-state solution, would provide the starting point 
for the next perturbation. This may well be true provided that at each step the steady motion is both 
statically and dynamically s .Die, and that the actual disturbances, e.g., the amplitude of oscillation, 
remain small. However, when the angle of attack exceeds a certain critical value o cr at which the steady 
motion is no longer stable, the linear theory predicts an exponential growth of the perturbation with time 
and, therefore, must cease to be valid after finite time. The motion of the aircraft under these conditions 
can only be studied by means of a nonlinear theory. 

In this paper we investigate three types of problems in which a steady motion becomes unstable. 

Padfleld (Ref. 2) studied a similar class of problems using the method of multiple scales, which is valid 
only for weakly nonlinear oscillations. We shall study the problem using bifurcation theory. This allows 
us to draw on recent mathematical developments (see Ref. 3) that are particularly well-suited to investigat- 
ing fundamental questions in linear and nonlinear stability theory. A numerical scheme based on bifurcation 
theory was proposed earlier (Ref. 4) for analyzing aircraft dynamic stability in a general framework. More 
recent work by Guicheteau (Ref. 5) demonstrates the considerable potential ?f bifurcation theory in flight 
dynamics studies, particularly toward establishing a method for the design of flight control systems to 
ensure protection against loss of control. While acknowledging the Importance of the aerodynamic model in 
determining aircraft stability characteristics, neither study contains an adequate assessment of the model 
requirements. The treatment of unsteady flow effects receives no attention. In contrast, we shall focus on 
just this aspect of the motion analysis. 

We shall restrict this study to the single-degree-of-freedom motions, e.g., pitching or rolling, of an 
aircraft about its trimmed flight condition. This enables us to analyze motions for which complete aerody- 
namic Information (exact analytical or numerical solutions) is available for certain shapes. The informa- 
tion is obtained from solutions of tho unsteady Inviscld flow equations (Refs. 6-12) or from results of 
experiments (Ref. 13). In this way it is possible to establish a form revealing a precise analytical rela- 
tionship between the basic aerodynamic confidents and the characteristics of the motion. 

Specifically we shall consider the following three examples t 

A. Pitching supersonio/hypersonic airfoil in rectilinear flow (Fig. 1); 



B, Flap oscillations In transonic flow (Pig. 2); 

C. Wing took of a slender delta wing In subsonic flow (Fig. 3). 


2. MATHEMATICAL FORMULATION 


Let the alroraft or flap be in level, trimmed flight until tlma t - 0 whan It la parturbad from Its 
trim position. During tha aubsaquant motion tha oantar of gravity oontlnuaa to follow a raotllinaar path at 
oonatant valooity V*. For a a Ingle-degree-of -freedom oaolllatory motion, tha aquation of motion ia 


d 2. 

I 3- 1 . Q( t;X) 
dt* 


Cl) 


whara 5(t) ia tha instantaneous angle-of-attack parturbation in axampla A, flap-angle perturbation In 
example 8, and roll angle in example C; I la tha appropriate moment of inertia; and G(t) la the corre- 
sponding Instantaneous perturbation moment of tha aarodynamlo forces. 

In Eq. (1), X repreaenta a sat of parameters defining steady flight at tha trim condition. Flight 
Maoh number M., ratio of specific heats V, and trim angle of attack g n are included in these param- 
eters. We shall consider X to be the trim angle of attack o m in examples A and C, and tha mean flap 
deflection angle in example B. In other words, all other defining parameters will be held fixed whan 
considering the consequence of varying X on aircraft motion characteristics. We assume that tha moment 
required to trim the aircraft or flap at X has bean accounted for, so that C(t;X) is a measure of the 
perturbation moment only. 

The instantaneous motion state, £(t) and £(t), and the instantaneous moment G(t) are a result of the 
interaction of aircraft or flap motion and the unsteady aarodynamlo forces from time zero to time t. 
Consequently, the instantaneous moment G(t) depends not only on tha Instantaneous motion state, g(t) and 
J(t), but also on the past motion history from time zero to t. This is to say that G(t) is a functional 
of £(tj)» (0 i t 1 S t) as described in Ref. 14. Thus 

t 

0<t) . ofeu,)] (2) 

0 

For motions for which C(t^) is analytic, the functional 1s equivalent to a function of an infinite set of 
variables, 

t 

G(t) • Gfttt,)] 

0 


let), ‘itt) 


<i n e(t) 

at" 


(3) 


For most problems encountered in the study of dynamic stability, tha motion i r slow although its amplitude 
may be finite or large. Under these conditions £(t), *C( t) , ... in Eq. (3) may be neglected and, as a 
first approximation 

G(t) - CU(t),£(t)) (4) 

We further assume that G ia r , analytic function of (and £. Expanding Eq. (4) as a Taylor series, 
we neglect terms of 0(£ 2 ) and b'. # ner for alow motions, and after reintroducing the parameter X, 


Q(tiX) - Oq((( t) ; A) 


iiUi o jttlthx) 


( 3 ) 


where t lsa characteristic length and 

G 0 (0»X) • 0 (6) 

as required at the trim condition. The form of the instantaneous moment G(tsA) in Eq. (8) Is consistent 
with the exact analytic solution for the pitching moment in example A (Ref. 10), the numerical solution Tor 
the hinge moment in example B (Ref. 12), and the experimentally derived empirical formula for the rolling 
moment In example C (Ref. 13). 
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In summary, the a lngle-degree-of -freedom motions considered in this papar will be based on the follow- 
ing mathematical problem 


2 • 

I * °o (| 5 ' x) * T °i (5|X) * 

dt • 

.6(0) - 6g 
6(0) - 6, 


(7.) 

(7b) 

(7o) 


The funotions G 0 UfX) and GjUsX) are generally nonlinear in i and have to be determined from the study 
of unsteady aerodynamics, either theoretically or experimentally. Evidently G 0 is related to the restor- 
ing moment and G 1 Is related to the damping moment. 

In many situations, it is known that when the parameter X reaches some critical value X qr , the 
aerodynamic damping G, vanishes and steady flight at the trio condition X Qr loses its stability. We use 
bifurcation theory to determine the motion characteristics of an aircraft or flap whose trim condition is 
near or beyond X Qr . 


3. BIFURCATION THEORY 

We Introduce the dimensionless time t • V a t/t, where 1 is a characteristic length, equal tc the 
chord length of the airfoil in example A, the ohord length of the forebody plus flap In example B, and the 
chord length of the wing in example C. Hereafter we use (*) to denote d/dt. We further let 

P(6.6sX) - G(6.6i*)/vl vf/t 2 ) 

• F 0 (6;X) ♦ 6F,(6iX) (8) 


Then Eq. (7a) may be written 




(9a) 


H • F(E.EiX) 


(9b) 


An expansion of F(S,5;X) in a Taylor series in £ and i and a ohange of notation u^ « S, u g • l 
yield for Eqs. (9) 


“I ■ V A)u j * B ijk a)u A * c ijki U) W» * 0(| “ |l,> • <l ■ 1,?) 


( 10 ) 


where 


A - 


[° '1 

L-S(X) “D(A)J 


(It.) 


SU) • -F£(O t A) , D(A) • -F 1 (OiX) 


(11b) 


V 0 * 


1 3 2 F 


2Jk " 21 3Uj 3u k 


(lie) 


u«0 


and 




1 3 J F 

c ijkt * 0 • C 2jki ’ JT a Uj au, K au t 


(lib) 


u-0 


(Although Eqs. (9) have been derived on the assumption of slow oscillations, subsequent bifurcation analysis 
of Ec. (9) will hold for general FU.EiX), i.e., as if no restriction had been placed on the magnitude 
of ?.) 


In Eqs. (ID, the tensors B and C represent the effects of ^nite amplitude to the second and third 
order. The following symmetry properties hold? 
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S 2Jk " S 2kJ 

C 2Jkt ‘ C 2J (k ' C 2k*j * C 2tjK ‘ C 21KJ 


On the basis of Eq. (to), 


wa discuss the llnaar and nonllnaar stability of the motion. 


(12a) 

(12b) 


Llnaar Stability Theory 


Tha a tab ill ty of ataady motion at tha trim oondltion 
the nature of tha alganvaluaa of tha matrix A. They are 


to Infinitesimal disturbances la determined by 


n,U) • j C-D(A) t /d 2 (a) - 4S(A) ] 


(13) 


unstaSr 1 ' SU) < In thlS 0aae ' * °’ n 2 < °* Th# «•»«* notion at thla oondltion A is always 
Case* li t S(a) > 0 . 


Ha : 0(A) < 0. 

Ilb t 0(A) > 0. 


In this case, Re(n t ) > 0 and the steady motion at A 
and the steady motion at A 


In this case, Re(n,) < 0 


is unstable, 
is stable. 


the ataady motion at tha trim cr.dltlc A inflllt iTJ?* darlvatlva D <»> Positive, is 

(Ref. 3) can b. used to show that stability of tc l " fin “” lnal disturbance. In fact, stability theory 
disturbance is sufficiently small. * * * dy motion in this case is assured if and only if the 

r™. g. - . »«■««.. ««. * « «. lu «, 

Sirs irrsr* ^ 

Thu ehihm.,, . 3 


ary A • a 


cr 


where 


studied presently 


SU cr ) > 0 


and 


DU or ) 


sbauiiiL; DuUnu 

The stability oharaotarlsttos near this boundary will be 


3.2 Hopf Bifurcation Theory 


At the dynamic stability boundary A - A cp , we have S(A ep ) > 0, hanca. 


W 


ti/snzr 

cr 


tlu. 


( 1 «) 


B»a exlatenca of purely imaginary eigenvalues of the matrix a at 1 . 1 <. . 

Hopf bifurcation (Refs. 3 and 15), slanalina a ehAn .. nu ._ , * at , X x er ls th * oharaotaristio sign of a 
crossing A . A cp . the st\£ «otlon mt h.d b..n f!.b^ to r w * **?* "° tl0, ‘ t0 P * ri0(llc 
resulting (after a transient motion has died away) in the e*iati,a.*25 btcoot * unatabla to disturbances, 
is periodic, in the vicinity of A • A tha oireii J ha ,xl8tenoa of • p » w notion, which (if it is stable) 

y er * n u«i r frnqueney of the periodic motion ls nearly equal 


For A slightly larger than 


In this soction we 


X or , the sigenvalues of the matrix A 


are 


1 


n, •-£ 0(A) t iQ(A ) 
2 


(15) 


where 


80) • /S(A) - D*(A)/t 

He shall assume that 


(16) 


D'(A ) < 0 
or 


le the usual case in applications. (The orse D'O ) > 0 
normalized eigenvector ;(A) associated with the eigenvalue 


(IT) 

can be treated in exactly the same way.) 
n(A) is 
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whereas the adjoint eigenvector c # (x) with eigenvalue ’n(X), which la f ihe oomplex conjugate of 



A. Hopf Bifurcation. The bifurcation solution u(t.X) nay be written as 

u • a(t)c ♦ a(t)c 


(18) 

n(x), is 

(19) 


(20) 


Following Iooaa and Joseph (Ref. 3. P» 125), we get 


a • eb 1 (a) ♦ e 2 b 2 (a) ♦ c 3 b 3 (a) ♦ 0(c 4 ) 
2 4 

s • [w Q ♦ c o» 2 ♦ 0(c )]t 
x * x er * *\ * 0<e4) 


( 21 ) 


(Ref. 11). The solution Is 


where, for brevity, we omit the lengthy solutions for e, b n , wj.^and X 2 
periodic in t with circular frequency equal to wq * e*«2 4 0(« )• 

B. Stability of the Bifurcation Periodic Solution. According to Floquet theory (Ref. 3)* the jtabll- 
lty of the bifurcation periodic solution Eq. (21) Is determined by the sign of an Index u. To 0(e ), u 
has the form 


D*(X U,e 2 ♦ 0(e ) 


( 22 ) 


The bifurcation periodic solution is stable if u < 0 and unstable if u > 0. Since we have assumed 
0*(X or ) < 0, stability depends on the sign of X 2 , with X 2 > 0 denoting stability and X 2 < 0 instabil- 
lty .It remains to cast X 2 in more recognizable terras. Considerable manipulation yields 


u 




2 

rF 

du 1 3 u 2 



(23) 


in terras of F(u 1P u 2t X). From Eq. (8) we see that the function F is directly related to the moment 
G(E,E;X) acting on the aircraft or flap which is performing a finite-amplitude oscillation £ around the 
trim position X. Equation (23) demonstrates that the stability of periodic motion near the dynamic stabll 
lty boundary X la determined by the behavior of the aerodynamic response 0(5. CiX) in that vicinity. 

e2 ’* 

With the assumption of slow oscillations under which the form of Eq. (5) was derived (terras 0(£ *0 
neglected), we may substitute Eq. (5) with Eq. (8) into Eq. (23) to get 


u 


7T < F 3 F i * Wi-o 




x-x 


or 


(24) 


Since from Eqs. (I 1 *) and (lib) 

-o • /5T V y • / ’ F 6 (0,x cr > 


Eq. (24) becomes, after using Eq. (8), 



(25) 


(26a) 


(26b) 


26-6 


Using Eq. (10) we get 



(26o) 


We have thus established the following orltoriont blfuroatlon perlodlo notion la stable or unatabla aooord- 
lng to 

u < 0 or u > 0 (27a) 

or alternatively, slnoa S(X op ) > 0, aooordlng to 

C 2112 (A cr ) < 0 «• C 2112 U or ) > 0 (27b) 

The two possibilities are well Illustrated In the fora of blfuroatlon diagrams as shown In figs. 4a and 4b. 

In a bifurcation diagram, the absolssa represents the blfuroatlon parameter 1, while the ordinate Is a 
oarameter characteristic of the blfuroatlon solution alone. In this oase It la c, a measure of the ampli- 
tude of the blfuroatlon periodic solution. Stable solutions are Indicated by solid lines and unstable 
solutions are indicated by dashed lines. Thus over the range of the blfuroatlon parameter X < X Qp where 
the steady-state motion Is stable, c Is sero, and the stable steady motion Is represented along the 
abscissa by a solid line. The steady motion becomes unstable for all values of x > x Qp as the dashed line 
along the abscissa indicates. Perlodlo solutions bifurcate from X « X 0P , either super or ltloally or 
subcrltlcally. 

When C 2112 (x cp ) < 0, hence u < 0 (implying X 2 > 0), the bifurcation Is called supercritical and its 
characteristic form is shown in fig. 4a. In this oase, stable perlodlo solutions (solid curves in Fig. 4a) 
exist for values of X > X Qp . The amplitude of the periodic solution at a given value of X - X 0P is 
proportional to c, hence Is vanishingly small when X - x gp is small, varying essentially as (X - X^) 1 ' 2 . 

When ^21 is^or) > 0 twice u > 0 (Implying X 2 < 0), the bifurcation Is called subcrltlcal and its 
characteristic form is shown In fig. 4b. In this oase, perlodlo solutions exist for values of X < X Qr , but 
they are unstable (dashed curve in Pig. 4b). The existence of stable periodic solutions for X > x Qr 
depends predominantly on the behavior of the damping Q^(CiA) for X > X Qp . If no stable perlodlo solutions 
exist for X > X QP , then when X is increased beyond X op the alroraft or flap may undergo an aperiodic 
motion whose departure from the steady motion at X « X QP Is potentially large. 

In the more likely event that stable periodic solutions do exist for X > X pp , their amplitudes must be 
finite, and not infinitesimally small, even for small positive values of X - x Qp , It Is likely that this 
branch of stable periodic solutions will join that of the unstable branch as illustrated in Fig. 4b. In 
this event, the form of the bifurcation curve for values of X < X op helps explain the situation where the 
steady-state motion could be stable to sufficiently small disturbances but become unstable to larger distur- 
bances. For x < x Qp , Fig. 4b suggests disturbances with small enough amplitudes (lying below curve 08) 
will die out and the steady motion will remain stable. However, disturbances with amplitudes sufficiently 
larger than those of the unstable branch may actually grow up to the ultimate motion aa t * •, which will 
be that of the stable branch of perlodlo solutions (curve BA in Fig. 4b). Finally, we note that if the 
motion does attain the stable branch of perlodlo solutions (say, for X < x pp ), then hysteresis effects will 
manifest themselves with further changes In X. When X Is increased beyond X QP , motion will continue to 
be periodic with finite amplitude (point A In Fig. 4b). If X Is now decreased below X Qpl periodic motion 
will persist, even at values of X where previously there had been steady motion when X was being 
Increased. Not until X is diminished beyond a certain point (point B in Fig. 4b) will the motion return 
to the steady-state condition (point C In Fig, 4b) that had been experienced when X was Increasing. 


4. APPLICATIONS TO AIRCRAFT MOTIONS 

In this aeotlon we shall apply the theory developed in Secs. 2 and 3 to three different aingle-degree- 
of -freedom motions. 


4.1 Pitching Motion of Supersonlc/Hypereonlo Airfoil 

We retain the symbol £ In all three oases to designate the perturbation variable. The other symbols 
win be replaced as necessary by the parameters relevant to the particular motion. 

In the oase of pitching osolllatlon (Fig. D of a aupersonie/hypersonlo airfoil in rectilinear flight, 
the parameter X la the trim angle of attack, designated in Fig. t, measured relative to the horizon- 
tal velocity vector. The instantaneous argle of attack la o(t), also measured relative to the horizontal 
velocity vector, so that the perturbation variable £(t) Is the inclination o(t) - o m from the fixed trim 
anile of attaok. The perturbation mernent Q in this oase Is the perturbation pitching moment about the 
center of gravity. Thus, 
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0< 1 1 A ) - qSlC^tio^.M^Y.h) 


(28) 


where* q La the dymunlo pressure, and S and l are the reference araa and length. The Inatantanaouo 
pitching-moment ooafrtoiant C m la alao a funotlon of tha trim angle of attack o m , the flight Mach num- 
bar M., tha ratio of apaotflo haat Y, and tha (dlraanalonlaaa) pivot axis poaltlon h. 

For large-amplitude alow oscillations, axact aolutlons for C m axlat for almpla ahapaa (Refs, 10 
ana n), and thair form la oonalatant with Eq. (5), 


C b ( ti V • c m 0 ( « (T),0 « ) * 4c « 1 (5(T,,a B > ' C m 0 (0 'V 


<29) 


Moreover, it la ahown that (Refs. 10 and 11) 


c . l (t,a « ) 


V* * a m> • 


(l ■ 0, 1) 


(30) 


Accordingly 


3C 3C_ 
“l "l 


a ■ 0, 1) 


(3D 


and Eq . (26b) reduces to 


' * ‘"bF^J . .0 

J m or 


(32) 


where the stiffness derivative S(o m ) and the damping- in-pitch derivative D(o ffl ) are related to the 
pitching-moment coefficient C m by 


S(oJ • -C* (o m ) 

m bIq (n 


D<0j 

m 


-C_ 


(33) 


and 


D(o > 
or 


(3U) 


A typical example of the variation of damping derivative D(o n ) veraua o ffl is shown in Fig.. 5 (Ref. 16). 


We conclude that when the trim angle of attack o. 


la increased beyond a critical value o or at which 


the aerodynamic damping vanishes, l.e., 0(o or ) ■ 0, the steady motion at o or loses its stability. A new 
finite-amplitude periodic motion, called the bifurcation solution, will replace it for values of o fc o Qr 
if it is stable. Stability of the bifurcation solution depends on whether D • ( c m ) /S( o m > is increasing or 


decreasing on crossing the stability boundary o or . 
By utilizing an approximate relation (Ref. 17) 


0( o ) 

m 


bCS(o ) 
cr 


S(a )] . 
m 


b > 0 


in Eq. (32), we get 


u 



o-o 

B 


or 


For the oaae of a flat-plate airfoil in supersonic/hypersonic flow, tha stiffness and damplng-in-pitch 
derivatives S(o n ) and D(o a ) are known exactly in analytical form (Ref. 7) for values of o ffl up to the 
shock detachment angle. It is shown in Table 1 that u > 0 for all combinations of M. and h. Thus, the 
bifurcation is subcrltlcal. The bifurcation perloulo solution la therefore unstable, which implies one of 
two alternatives. The steady motion is either replaced by a finite-amplitude periodic motion accompanied by 
hysteresis phenomena, or by a potentially large aperiodic motion. In either case the loss of stability of 
the steady-state motion must be accompanied by a discrete change to a new stable state. 


M .2 Flap Oscillation in Transonic Flow 

The case of a flap oscillating about a hinge is conceptually similar to that of an airfoil pitching 
about a pivot axis. In effect, the oscillating flap may be regarded as an airfoil pitching in the nonuni- 
form onooming stream that rssults from s uniform free stream passing over the fixed forebody (Fig. 2). 
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In this case (Fig, 2), <5( t ) 1a the instantaneous flap deflection angle, measured relative to the under 
flaotad position. Tha parameter X here la tha flxad main deflection angXa « m , aXao maaaurad from the 
undaflaotad position. Thus, $(t), tha perturbation variable, is tha InoXination 6(t) - fl m from tha flxad 
naan daflaotlon angla <$ m . Tha parturbatlon momant 0 Is tha parturbatlon hlnga moment about tha hinge 
line. Thus 

0(t»X) - qSlC^( t (35) 

where C h is the instantaneous hinge-moment coefficient. Xn direct parallax with the oaolllatlng airfoil 
case of Eq. (29), we assign C h the form 

• y^'V * < y ,<tM - > ‘ V 0 ’ 4 -’ <36) 


Equation (36) is consistent with Eq. (5). A similar form was validated in Ref. 12 by comparison with 
results of large-jcale numerical integrations of the coupled Inertlal/flow field equations. 

Equation (5) or Eq. (36) for the hinge-moment coefficient allows calculation of C h ^ (g( x) \6 m ) by 

linearizing the flow field equations about the steady flow corresponding to fixing the instantaneous deflec- 
tion angle, S(t) ♦ d m (see Eq. (39)). Calculation of (^( t) 1 6 m ) from the linearized equations may 

require less computing time than the nonlinear method used in Ref. 12, yet it Is consistent with the level 
of approximation leading to Eq. (5) or Eq. (36). This simplification of the calculation of < 5( x) ; d m ) 

reinforces the main conclusion of Ref. 12 that the coat of computations to determine the flap's motion will 
be very small if the aerodynamic contribution to the equations of motion is modeled, as compared to the 
computational cost required to solve the coupled inertlal/flow field equations. 

It was further established (Ref. 11) that 


c h ■ c hi <« ♦ V • (i • °’ 15 


(37) 


which in turn reduces Eq. (26b) to 


where 

SUJ ■ -Ci <«_) , 0(6 ) • -C. (« ) , ) • 0 (39) 

m h Q m m m or 

We therefore reach a conclusion regarding the characteristics of the flap oscillation near the mean flap 
deflection angle d op whioh is direotly analogous to that of Sec. 4.1 regarding the characteristics of the 
airfoil pitching motion in the neighborhood of o Qr . 



Numerical results from Ref. 12 for S(6 n ) and D(8 B ) are reproduced (by means of spline fitting) in 
Fig. 6. These are results for an NACA 64A010 airfoil at zero mean angle of attack in a transonic stream 
with • 0.8 and Y » 1.4. It is found that 

[ a D *<6 _>1 

^^ vv ”"* 20 “ <cr ’ 3 ' 7 ‘ 

Henoe, the bifurcation is subcrltlcal and the bifurcation periodic solution is unstable. Just as we found 
for the flat-plate airfoil in supersonic/hypersonic flow, loss of stability of the steady-state transonic 
flow about the flap at 6 m • 6 0r implies a discrete change to a new stable state. The new state may be 
either a finite-amplitude periodic oscillation about 6 n accompanied by hysteresis phenomena, or it may be 
a potentially large aperiodic departure from the steady state. The fact that D(6 n ) becomes positive again 
at |*J • 17° (Fig. 6(b)) implies that the first of these alternatives will be the preferred one (Ref. ID. 


4.3. Wing Rock of Slender Delta Wing in Subsonic Flow 

The phenomenon known as "wing rock" of a slender delta wing in subsonic flow has been the subject cf 
Intensive Investigation by many researchers (Refs. 13, 18-21). It is now well-documented (Refs. 13 and 18) 
that as the trim angle of attack o m is raised incrementally, a critical value o cr is reached where roll- 
damping moment changes from stabilizing to destabilizing. Ericsson (Ref. 2D has shown how the leeside 
vortioes, which appear and grow stronger as o m is raised, may lag behind a rolling motion of the wing. 

The pressures they Induce may contribute a rolling-moment component proportional to rolling velocity that is 
destabilizing and eventually becomes large enough to surpass the stabilizing component that had existed 
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alorw at mnallar a m . Consequently, at o m £ a , small disturbances in the flow field exalte a growing 
roll oscillation. Bifurcation theory again can be applied to help determine the ultimate state of the 
oaoillation. 

For this case, we take the parameter X to be the fixed trim angle of attack o m » while the perturba- 
tion variable C(t) is the roll defleotlon angle ♦(t) (Fig. 3)» measured from the position of aero roll 

angle. The perturbation moment 0 la the rolling moment about the wing longitudinal axis. Thun 

G( t | X) - qStC ft (t*o m ,A) (40) 

where, for incompressible flow, the instantaneous rolling-moment ooeffloient depends on the trim angle 
of attack o m and the aweepback angle A of the delta wine. Results of experiments (Refs. 13 and 18) and 

numerical computations (Ref. 19) have shown that, for a given — A, it is a good approximation to the instan- 

taneous rolling-moment coefficient to take 

c i - v e(t),, '» ) + ic » 1 U(t),0 » ) (,,1) 

which Is oonsistent with Cq. (5). 

As an example, for the 80° sweep-baok flat delta wing studied in Ref. 18, is approximated by a 
power series (Ref. 19) which leads to the following equation for the rolling motion 

i « F U4io m ) . 

■ [yo m H ♦ b 3 (o 0 >5 3 ] ♦ to 0 ♦ b 2 Co m > ♦ b 4 (o ra )5 2 3 

' F o (e,# l n ) + ( " 2) 

Taking account of the scaling factors ^adopted in Ref. 19, 

W * K2c i a i (o m > * ° ‘ ’• 3) (43a) 

bj(° m ) - <C 1 a J (o m ) , (J * 2, 4) (43b) 

In Bqs. (43), k is a factor accounting for the difference between scalings adopted for the time variable in 
Ref. 19 and this study. From Ref. 19, k - 2c/L - 2 x 0.429/0.107 » 8, C 1 • 0.088. The a*(o m ) are tabulated 
in Table 5 of Ref. 19, which yields the following table for the b i (o m )i 


o 

3 

O 

o 

15® 

o 

O 

CM 

25° 

b, -0.0265 

-0.0721 

-0.1977 

-0.3320 

b r -0.01C1 

0.0090 

0.0596 

0.0959 

bj -0.1222 

-0.2714 

-0.0501 

0.2894 

b 4 +0.1491 

0.1159 

-0.1799 

-0.9977 


From Eq. (11b), the linear contribution to the damping in roll is 


- -V°'V - - b o- w 


(44) 


where -b Q > 0 la proportional to the wind-off roll-damping moment due to bearing friction in the experi- 
mental setup used in Ref. 16. In order to compare our theoretical prediction with the experimental results 
of Ref. 18, we choose a value of -b 0 such that Eq. (44) yields zero damping in roll at the same angle of 
attaok, o or - 18.6°, as that reported in Ref, 18. The funotion D(o m ) Is then plotted in Fig. 7 by means of 
spline fitting. 


From Eq, (26o), the index for stability of the bifurcation periodic rolling oscillation is 
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(45) 


At « cr • 18.6*, we find from the table above that bj(o or ) - -0,1591 < 0. From the spline-fit curve for 
bj|(o cp ) (Fig. 8), we find that b^(o Qr ) • -0,05473* Hence u < 0 and the bifurcation is supercritical, 
implying that the bifurcation periodlo roll oscillation is stable. This is in agreement with the experimen- 
tal finding in Ref. 18. 
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W« now further oomp.ro th. superartttoal bifurcation diagram predicted by th* prM.nt theory with 

WU-nux result a from Ref. 18. Combining Cq. (21o) with Eq. (22) to .llmln.t. i 2 and thw/ualng' 
&q . ( HD) , W© get * 


Q m “* " "C 

m • op 


* r-« “< i 

b «rj v 

H., 
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(46) 


Tho laat^aquatlon was obtaln.d by u.lng Eqa. (25) and (42). From Fig. 7 (noting that 1" - t/180), w. g.t 
d (o or ; ■ * " CN 


"0.61 31* Hence Eq. (46) reduoes to 


°m 


or 


0.1120 c 


(47) 


whera c is the amplitude of tho bifurcation periodic aolution (Ref. 11). Equation (47) la plotted In 
* £? r C !2 e 9weap '' back dalta wing and compared with the experimental reaulta of Levin and Katz 
(Ref. 18). The agreement la excellent near o or , where the bifurcation theory la particularly applicable. 

It should be noted that although u < 0 and the bifurcation ia supercritical , the magnitude of |u| 
la small in the example. This Implies that the bifurcation la close to the boundary between aubcritical and 
supercritical. Consequently, on Increasing the angle of attack o n by a small amount, the amplitude of the 
resulting periodic motion will be quite large. This ia confirmed by the experimental results ahown in 

with an amplir^ Pl f , 22 n 6« eaSln8 ° m ^ P * at ° or ( * l8,60) re8ults a stable periodic roll oscillation 


5. CONCLUDING REMARKS 

We have shown how bifurcation theory can be used to study the nonlinear dynamic stability characteris- 
tics of an aircraft or flap subject to aingle-degree-of-freedom motion about its trim position. When the 
bifurcation parameter X (e.g., the angle of attack) is increased past the stability boundary X- r , where 
the aerodynamic damping vanishes, the steady motion loses its stability. This results in a finite-amplitude 
periodio motion after the transient motion has died away. We have also established a simple criterion for 
the stability of the bifurcation periodic motion in terms of the aerodynamic coefficients. The theory 
predicts that bifurcation solutions are unstable (aubcritical) for the pitching airfoil in 
supersonio/hypersonic flow and for flap oscillations in transonic flow. Bifurcation solutions are stable 
(supercri tical ) for roll oscillations of the slender delta wing in subsonic flow. The latter prediction is 
in good agreement with available experimental results. 

When the theory predicts subcritical bifurcation, the abrupt change in the motion that results when 
X is increased past X cr may cause an abrupt structural change of the flow field. This in turn may render 
invalid the form of the perturbation moment of the aerodynamic forces that was used. Under these condi- 
tions, aerodynamic information in a different form may be required. However, the theory as developed in 
this paper is valid up to X cr and can be used to predict the onset of subcritical bifurcation. 

When the theory predicts supercritical bifurcation, the bifurcation periodic solution is stable to 
small disturbances for a range of X beyond X op . With further increase in X, however, the periodic 
motion may lose its stability. This in turn may result in another bifurcation at X • A* > x rtnl which could 
be either subcritical or supercritical. The reaulting bifurcation solution may be quasi-periodic, and the 
sequence of bifurcations may continue. 

Finally, the theory may be generalized to apply to aircraft motions involving more than a single degree 
of freedom. For motions involving two degrees of freedom or more, one should be aware of thu possibility of 
chaotic behavior occurring after a finite number of successive bifurcations. 
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Fig. 4. Typical forms of bifurcation diagrams 
near the dynamic stability boundary X where 
DU cr ) • 0. (a) Supercritical, u < 0. (b) 

Subcritical, u > 0. 



Pig. 5. Damping-in-pltch derivative 0 versus 
angle of attack o a for a biconvex circular arc 
airfoil (Ref. 16), ft • 0.5, Y ■ 1.4, t • 0.075. 



Fig. ba. Stiffness derivative S versus mean 


deflection angle « m of a 
64A010 airfoil (Ref. 12). 

0m * 0. 


transonic flap on NASA 
M. - 0.8, Y • 1.4, 



Fig. 6b. Damping derivative D versus mean 
deflection angle 6 m of a transonic flap on NASA 
64A010 airfoil (Ref. 12). M. - 0.8, Y - 1 4 
°m • «or " 2 *? # * 



Fig. 7. Damping in roll D of an 80° sweep-back 
flat delta wing versus angle of attack o m in 
incompressible flow (Refs. 18, 19), Eq. (44). 

°cr » 18.6*. 
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